function d_filt = filterBank(d)
bank = [[0.001 4];[4 8];[8 12];[12 15];[15 30];[30 95]];
energy = zeros(size(bank,1),1);
d_tot = d*0;
for f=bank'
 d_filt = myFilter(d,f,1);
 d_tot = d_tot+d_filt;
 energy(find(ismember(bank(:,1:2),f','rows') == 1)) = sum(d_filt.^2);
end
[v,I] = max(energy);
d_filt = myFilter(d,bank(I,:));
plot(d,'b');hold on;plot(d_tot,'r');
end